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ABSTRACT 

A comprehensive new approach is presented for deriving probability densities of physical 
properties characterizing the lens and source that constitute an observed galactic microlensing 
event. While previously encountered problems are overcome, constraints from event anoma- 
lies and model parameter uncertainties can be incorporated into the estimates. Probability 
densities for given events need to be carefully distinguished from the statistical distribution 
of the same parameters among the underlying population from which the actual lenses and 
sources are drawn. Using given model distributions of the mass spectrum, the mass density, 
and the velocity distribution of Galactic disk and bulge constituents, probability densities of 
lens mass, distance, and the effective lens-source velocities are derived, where the effect on 
the distribution that arises from additional observations of annual parallax or finite-source ef- 
fects, or the absence of significant effects, is shown. The presented formalism can also be used 
to calculate probabilities for the lens to belong to one or another population and to estimate 
parameters that characterize anomalies. Finally, it is shown how detection efficiency maps for 
binary-lens companions in the physical parameters companion mass and orbital semi-major 
axis arise from values determined for the mass ratio and dimensionless projected separation 
parameter, including the deprojection of the orbital motion for elliptical orbits. Compared to 
the naive estimate based on 'typical values', the detection efficiency for low-mass compan- 
ions is increased by mixing in higher detection efficiencies for smaller mass ratios (i.e smaller 
masses of the primary). 

Key words: gravitational lensing - methods: statistical - binaries: general - planetary sys- 
tems - Galaxy: stellar content. 



1 INTRODUCTION 

During the recent years, more than 2000 microlensing events have 
been observed and corresponding model parameters have been pub- 
lished. However, these model parameters in general do not coincide 
with the underlying physical characteristics of lens and source star, 
which are their distances from the observer, Z?l and Ds, respec- 
tively, the mass M of the lens, and the relative proper motion /^ls 
between lens and source. For 'ordinary' events, compatible with 
rectilinear motion between point-like sources and lenses, the only 
parameter related to these characteristics is the event time-scale 
tE = Oe/^ls, which corresponds to the time in which the source 
moves by the angular Einstein radius 

„ I4GM D s - D L UGM ttls 

relative to the lens, where 7r L s = 1 AU (Z)j~ 1 — Dg 1 ) denotes the 
relative lens-source parallax. 
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With the physical lens characteristics being statistically dis- 
tributed according to the mass density and velocity distribution 
of lenses and sources as well as the mass spectrum of the lenses, 
the distribution of observed parameters in the ensemble of galac- 
tic microlensing events can be use d to measure these distributions. 
iDe Ruiula. Jetzer & Massd ll99lh have shown explicitly how sta- 
tistical moments of the observed time-scale distributions translate 
into moments of the underlying mass spectrum of the lenses. 

A different question is posed by asking for the stochastical dis- 
tribution of physical lens and source properties given the observed 
model parameters for a single realized event. In the literature, this 
distinction has frequently not been made strictly enough, leading 
to some confusion. In particular, the probability density of the lens 
mass averaged over all observed events does not converge to the un- 
derlying mass spectrum. By quoting a probability for the lens mass 
in a given event to assume a specific value. lDe Ruiula et alj ll99ll) 
did not produce a meaningful result, since the probability for any 
random variable to assume a specific value is zero. Further com- 
mon misconceptions exist around a 'relative probability', which is 
not defined, and a 'most-probable value', which does not exist ei- 
ther. A finite probability can only be attached to a finite interval of 
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values, given as the integral over the probability density of the con- 
sidered quantity. iDominikl ll998al) realized that in order to derive 
information about the lens mass and other properties for a given 
event, one is dealing with a probability density, which carries the in- 
verse dimension of the quantity it refers to, rather than a likelihood 
(e.g. lAlcock et alJll995t> . In fact, likelihood functions and proba- 
bility densities are different entities, which can be seen explicitly 
from the following property: If a likelihood for a quantity a has a 
maximum at some value ao, the likelihood for any function f(a) 
of the quantity a has a maximum at /(ao), whereas such a prop- 
erty does not hold for probability densities, i.e. {/(ffl)) 7^ /((«)) 
ma y occur, where (a) deno tes the e xpectation value o f a. However, 
like lDe Ruiula et alj Jl99lh before. IPorninTkl Jl998di still failed to 
realize that a statistical mass spectrum of the lens population needs 
to be assumed along with the space- and velocity-distributions of 
lenses and sources. 1 Moreover, confusions around the ill assump- 
tion of a fixed-mass spectrum oc 8(M — Mfi xe d) led to incorrect re- 
sults for power-law mass spectra, where the power-law index would 
have to be shifted by unity in order to obtain correct expressions. 
For determining the event characteristics of MACHO 1997-BLG- 
41. lAlbrow etaflfeOOOiJ) used a discretization of the statistical dis- 
tributions of the basic lens and source properties in the form of 
a Monte-Carlo s imulation, which is a variant of the approach of 
|Pominiki ri998a). Some of the related ideas have been further de- 
veloped into a related formalism arguing on the basis of Bayesian 
statistics used in the analysis of OGLE 2003-BLG-423 <Yoo et alJ 
l2004l) . where the Galaxy model is used as the prior for the model 
parameters. 

In this paper, a revised comprehensive framework is presented 
for combining the model parameters as determined from the ob- 
servations with Galaxy models in order to estimate physical lens 
and source properties for a given event. This refined approach over- 
comes the previously encountered problems and allows the inclu- 
sion of model constraints from event anomalies as well as model 
parameter uncertainties. Moreover, by considering different lens 
populations, a probability that the observed event with its param- 
eters results from one or the other is obtained, which is taken into 
account for deriving the probability densities of the event charac- 
teristics. Rather than having to rely on Monte-Carlo simulations, 
all results are obtained in the form of closed expressions by means 
of integrals over the statistical distributions of the lens and source 
properties. 

In Sect. |2| the role of the event rate for deriving the desired 
probability densities of lens and source properties is discussed, 
while Sect.|5|looks at the global properties of the ensemble of mi- 
crolensing events such as the distribution of the event time-scale 
and the contribution arising from different lens populations. The 
probability densities of key properties of the lens, namely its mass, 
distance, and relative velocity with regard to a source at rest, that 
follow from a measurement of the event time-scale and the Galaxy 
model are derived and discussed in Sect. |4| whereas Sect. [5] fo- 
cusses on the implications if further constraints arise either from 
the measurement or from upper limits on additional model param- 
eters, where the two cases of annual parallax and finite source size 
are considered explicitly in detail. Sect.|6|then discusses probability 
densities of further quantities such as the parallax and finite-source 
parameters as well as the orbital semi-major axis and orbital period 
for binary lenses, before Sect.|7|shows how the presented approach 
can be used to determine the detection efficiency for companions 



By neglecting this, an implicit assumption is made. 



(such as planets) to the lens star as function of the physical proper- 
ties of the system. Sect, [^finally provides a summary. The under- 
lying probabilistic approach is presented in Appendix^] details of 
the adopted Galaxy model can be found in Appendix [B] and Ap- 
pendix|C|discusses the statistics of the orbits of binary systems and 
in particular the projection factor between the actual angular sepa- 
ration and the semi-major axis. 



2 PROPERTIES DETERMINING MICROLENSING 
EVENTS 

Microlensing relies on the chance alignment of observed source 
stars with intervening massive objects acting as lenses, where the 
degree of alignment is characterized by the angular Einstein radius 
as defined by Eq. Q, which depends on the lens mass M as well 
as on the source distance Ds and the lens distance Dl- With a 
two-dimensional angular separation 6 between lens and source, the 
magnification of the source star caused by the gravitational field of 
the lens in general depends only on u = 0/8e, while for a point 
source it even depends on its absolute value u — \u\ only, taking 
the analytical form (e.g. iPaczvnskl 19*81 



A{u) 



u -\- 2 
u y/u 2 + 4 



(2) 



The basic properties of point-like lenses and sources that af- 
fect the microlensing light curve are the source magnitude ms, the 
source distance Ds, the lens mass M, the lens distance Dl, the rel- 
ative proper motion between lens and source uls, taking into ac- 
count the motion of the observer, and the blend magnitude tub- If 
one considers the source distance Ds as well-constrained (there is 
no problem of including an uncertainty on this parameter as well), 
and uses the fact that there is no correlation between lens proper- 
ties and the source or the blend magnitude, it is sufficient to con- 
sider lens mass M, lens distance Dl and proper motion /ils as 
the descriptive properties for a given microlensing event. A binary 
lens involves further characteristics, namely its mass ratio q and 6 
orbital elements that can be chosen as the semi-major axis a, the 
eccentricity e, three parameters describing the orientation of the 
orbit (such as the inclination, the longitude of the ascending node, 
and the argument of perihelion), and finally an orbital phase (such 
as the mean anomaly at a given epoch). Distributions of the mass 
ratio q, the semi-major axis a, and the eccentricity e are pairwise 
correlated and also depend on the total mass M of the system and 
the actual types of stars involved, where our current knowledge on 
these is rather limited. 

Let v — Dl /^ls denote the effective velocity at the lens dis- 
tance that corresponds to the proper motion /ils, while the Einstein 
radius te = Dl 8e is the physical size of the angular Einstein ra- 
dius #e at this distance. With the mass spectrum $m (M) and the 
effective transverse velocity v being distributed as $„ (v), the con- 
tribution to the event rate by lenses in an infinitely thin sheet at 
distance Dl with masses in the range [M, M + dM] and velocities 
in the range [?;,« + dv] is given by 



dr = wo &m(M) v $ v (v) r E dM dv dD L . 



(3) 



where p(Dl) is the volume mass density, so that the differential 
area number density reads 



An = $ A f (M) dM dD L , 



(4) 
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and wo is a dimensionless factor representing a characteristic width 
that defines the range of impact parameters for which a microlens- 
ing event is considered to occur. Commonly, an 'event' is defined 
to take place if the source happens to be magnified by more than 
an adopted threshold value At, i.e. A At, where the choice 
At = 3/V5 w 1.34 corresponds to ht = 1, according to Eq. J3, 
which means that the source passes within the angular Einstein ra- 
dius of the lens, and therefore wo — 2. 

Instead of Dl, let us use the dimensionless fractional distance 
x = Dl/Ds, which is distributed as & x (x) = Ds p(xDs) /E, 
with E = J^ s p(D-l) cIDl being the total surface mass den- 
sity. Let us further assume that the mass spectrum $a/(M) is not 
spatially-dependent and involves a minimal mass M m in and a max- 
imal mass Mmax> while the velocity distribution & v (v, x) depends 
on the lens (and source) distance. With these definitions and as- 
sumptions, the event rate reads 




dM x 



v & v (v, x) dv y/x(l — x) <&x{x) dx 



(5) 



As discussed in detail in Appendix [A] the corresponding weight 
function 

Q(M,v,x) 



wq — vr E (D s , M,x) 



= Wo 



AG 1/2 E i— — 



(6) 



for the basic system properties a = (M, v, x) provides probability 
densities for any lens property that can be expressed as function of 
the basic properties by means of Bayes' theorem. 

By introducing a dimensionless velocity parameter £ = v/v c , 
where v c denotes a characteristic velocity, and with 
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(7) 



being the Einstein radius of a solar-mass lens located half-way be- 
tween observer and source (x = 0.5), the event rate can be written 

as 



' M max /M 



© 



\M min /M f7 



© 



M/Mq (M/Mq) 



d(M/M ) x 



C $ c (C, as) dC y/x(l-x) $ x (x) dx 



(8) 



with To = 2 wo r E q v c E and the dimensionless distributions 
*c(C) = v c $v{vcO and $ m/Mq (M/Mq) = Mq#m(M). 
With the definition of To, the weight function takes the form 



Q(M, C, x) = To (M/Mq) 1/2 C y/x(l-x) 



(9) 



mass lens and uniform motion of the lens relative to the line-of- 
sight from the observer to the source, the only observable that is 
related to the physical parameters of the system is the time-scale 



= 8e/pls 



AGM 



D s x{l 



(10) 



which thus involves all the basic properties M, v, and x. 
For an obtained best-fit estimate , let us define 



(0) 
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where r E q is given by Eq. Q. With depending on the basic 



system properties as 77^ = 2 JM/Mq y/xjl - x)/(, Eq. iA5l 
applied to the expression for the event rate T as given by Eq. jSJ 
yields the corresponding event rate density as 



/ (o)s _ 4 Tp 
"fvt E (Vt E ) - r (o) 



Af min /M0 
X <&M/Mq (M/Mq) X 



^M/Mq X 



x j $ c (2yM/M0 V /i(l-x)/r,|» ) , I ) x 



x [a;(l-x)] 3/2 $ ;c (s)da;d(M/M0). 



(12) 



while 7t E (t B ) = («c/r Ei 0)7 TOB [(t E u c )/r Bi 0] is the corre- 
sponding density of fE, and 7t E = 7t E /r gives the distribution 
of event time-scales arising from the lens population. 2 

If the lens may belong to one or another population with dif- 
ferent mass spectra, mass densities, and velocity distributions, the 
event rate density 7t E (£ E ) for the observed event time-scale 
provides a means to decide to which population the lens objects be- 
longs. Namely, the probability for the lens to be drawn from each of 
the populations is proportional to the corresponding event rate den- 
sity. Since the event rate density is proportional to the surface mass 
density along the line-of-sight, conclusions about the latter can be 
drawn, e.g. a likelihood for a certain surface mass density can be 
obtained on the assumption that the lens in the considered event 
with time-scale (or possible additional observables) belongs 
to a chosen population. Such considerations are of special interest 
with regard to the mass content of the Galactic halo and the still 
open question what fraction of the observed microlensing events in 
the direction of the Magellanic Clouds is cause d by lenses in the 
Magellanic Clouds themse lves (e.g lSahu &* Sahu 1998: lGvuk et alJ 
l200(j|Mancini et all2004l) . 

For a source located in the Galactic bulge, namely in the direc- 
tion of Baade's window with (l,b) = (1°,— 3.9°) at a distance of 
Ds =8.5 kpc and the lens residing in the Galactic disk or bulge, 
Fig. Q shows the distribution of time-scales and lens masses in- 
cluding the contributions of the individual lens populations, while 
Fig.|2|shows the fractional contributions Kt E (te) of disk or bulge 
lenses as function of the observed event time-scale t 



(0) 



disk 



7 t d r/(7 t d r+7t ulgc )^ 



bulge 



ft 



where 

/(7 t d r+7 t b : lse ), 



bul S c //..disk 



3 DISTRIBUTION OF PROPERTIES FOR THE 
ENSEMBLE OF EVENTS 

For ordinary microlensing light curves, i.e. those that can be ap- 
proximated by lensing of a point-like source star by a single point- 



2 Instead of eliminating the integration over by means of the 5-function, 
one can alternatively eliminate the integration over da; or d(M / Mq), 
which results in expressions that correspond to just two of the infinitely 
many possibilities to transform the remaining integration variables. 
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Figure 1. Distribution of event time-scales 4e and lens masses among all 
created microlensing events for a source located in the Galactic bulge at 
D$ = 8.5 kpc in the direction of Baade's window (l,b) = (1°,— 3.9°) 
with the individual contributions of disk and bulge lenses. Details of the 
adopted Galaxy model can be found in AppendixlBl 



and 7i E (£e) is given by Eq. IA5\ . TableQlists the fractional contri- 
butions K?' sk and K? ulgc for selected time-scales. All details of the 

C E tE 

assumed mass spectra, mass densities, and velocity distributions for 
the underlying populations can be found in AppendixlBl Among all 
created events, 35 per cent are caused by lenses in the Galactic disk 
and 65 per cent by lenses in the Galactic bulge. Only for timescales 
tE < 2 d, disk lenses provide a significantly larger contribution 
than bulge lenses, whereas the latter dominate for 2 d < £e J$ 40 d 
and for t B > 100 d. For 40 d < t E < 100 d, both popula- 
tions yield comparable contributions. One finds a median time- 
scale tE ~ 18 d, or tE ~ 17 d for bulge and 4e ~ 24 d for 
disk lenses. The distribution of jt E does not properly reflect that of 
the time-scales observed by the experiments, since their sensitivity 
for detecting an event depends on the event duration. In particular, a 
significant fraction of events with short time-scales is missed with a 
roughly daily sampling. The median mass is 0.32 Mq, with about 
1/3 of the events caused by lenses more massive than 0.5 Mq, and 
about 1/5 by lenses heavier than 0.8 Mq. 



4 PROBABILITY DENSITIES OF SYSTEM PROPERTIES 
FOR ORDINARY EVENTS 

4.1 Lens mass 

Let us define /io = M/Mo, where the characteristic mass Mo is 
assumed for the velocity v c and the lens being located half-way be- 
tween observer and source (x — 0.5). Hence, with r)^ as defined 
byEq. GQ,Mo = [vl°J] 2 M & , while nf n = [tiT^M^/Mq 

and M r x = [^'r'M^/MQ. 

With /io being related to the basic properties as /io = 
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Figure 2. Fractional contributions re? lsk and kJ 3 "'^ to the event rate den- 
's tE *E 

sity 7t E (4e ) as a function of the time-scale t-g of an observed event. For a 
few specific tE, the resulting values of re^ E sk and K^ lge are listed in Ta- 
ble0 The source is located in the Galactic bulge at Ds =8.5 kpc in the 
direction of Baade's window (I, b) = (1°, —3.9°). 



Table 1. Fractional contributions K^ lsk and re! 3 "' 60 of disk or bulge lenses 

*e *e fc 

to the event rate density 7t E (tE ) of a Galactic bulge source at Dg = 

8.5 kpc in the direction of Baade's window by lenses in the Galactic disk 

or bulge, respectively, for selected event time-scales tE ■ 



t E /(ld) 




bulge 


5 


0.25 


0.75 


10 


0.24 


0.76 


20 


0.35 


0.65 


40 


0.50 


0.50 


80 


0.52 


0.48 



< k = 7 t d E sk /(7 t d E sk +7 tE lEC ) and," = ^/ht^+y^l 
where 7t E is denned by Eq. <A5t and can be calculated by means of 
Eq. {12|. 
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2 (M/Mq), one easily finds with Eqs. jATBl and the 



probability density of no for an event with measured to be 



p (0) 
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4 Tn 



1m 



MO $ 



M/M 
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(0) 



$ c (^2 y/n x{l-x),a^ [x{l-x)\ i/2 $ x {x)Ax. (13) 



By means of Eq. <A16t . the distribution for a fuzzy value 
of then follows with p M0 (fj,o; Vt E ^) being evaluated for every 
corresponding Tj^ . Frequently, it is more adequate to represent 
the lens mass on a logarithmic scale. The probability density of 
lg (M /Mq) simply follows as 



Plg(M/M m ) (lg(M/M Q ),^ 



(OK 



(0) 



1 -2 M 

T7~ ln 10 Pmo 



M, 



© 



It 



(0) 



(14) 



2 M_ 

8.5 kpc in the direction of 



For a source located at Ds 
Baade's window with (l,b) = (1°,— 3.9°) as used in the pre- 
vious section and the lens residing in the Galactic disk or bulge 
with the Galaxy model described in AppendixlBl Fig. [3] shows the 
mass probability density Pi s (m/Mq) of \g(M /Mq) for selected 



values of the observed event time-scale t^ ' . In addition to the re 
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Figure 3. Probability density of the lens mass M for selected observed 
event time-scales tj;. With the event rate density 74 E (4e). the relative con- 
tributions by disk or bulge lenses are K d lsk — ~> dlsk / 1 ~. dlsk -l. ^, bul s c 



bulge / / disk 



7,t k /(7 t t k + % E n 

7 ( Jf sc ), respectively. The figures show 



"* '"t E - T tB ° /(7tE 1 n E 
the contributions Kt E Pi s (m/Mq) °f eacn population along with the total 

probability density Pi s (m/Mq) °f 1s(M/.Mq). For the chosen values of 

tE. me coiresponding fractional contributions Kt E (4e) are listed in TablelTI 



for M > 1 Mq forces the expected mass to be more narrowly 
distributed with tE rather than to follow the naive M oc t% law. 
In particular, the mass M — exp l0 [(lg(M / Mq)^] Mq spans 
only 1.5 decades between 0.09 Mq and 3 Mq for time-scales 
3 d ^ tE ^ 150 d, where the inclusion of one standard deviation 
extends this range to 0.03 Mq ... 15 Mq. 



4.2 Lens distance 

Similarly to the treatment of the lens mass, one finds the probability 
density of the fractional lens distance x for an event with observed 
4 0) (and related ?7 t (0) ) with Eqs. fA15l and dj as 



p! 0) (*;<*) 



(OK 



4 To 



'It 



(0) 



/Jm)®M/Mq (fJ-o 
X $ c ^2 y/ flo x(l - x), X^j dflQ 



x(l-x)] 3/2 $ x (x) X 



(15) 



4 r 



/Mq) 



lg(Mir 



(M/M ) 3/2 $ m/M0 (M/Mq) x 

lg(A/ min /M ) 

x 3> c (2 v^/M© Vx(1 - z)/^ 5 , 1) d[lg(M/M )] . (16) 

Fig. [4] shows the probability density p x of the fractional lens 
distance for selected time-scales, while the expectation value (a;) 
and the standard deviation o x for different are shown in the 
middle panel of Fig.|5|as well as in Table|2| As before, the source 
is assumed to be located in the Galactic bulge at Ds = 8.5 kpc in 
the direction of Baade's window, and the Galaxy model described 
in Appendix |B|has been adopted. Shorter time-scales favour larger 
fractional lens distances, while longer time-scales prefer the lenses 
to be closer to the observer, in accordance with the disk population 
yielding the slightly larger contribution to the event rate density 
for 40 d < 4e ;$ 100 d, whereas the disk dominates for smaller 
time-scales unless tE < 2 d. 



suiting P\ s (m/Mq) from both possible lens populations, their in- 
dividual contributions K tE P\^ 

are shown, where the factors are determined by the event rate den- 
sity 7tE (t E ) as K d E sk - - disk "- disk 1 ~ bulgc ^ ~ J - bulgo 



j bulge bulge 
(M/A/q) ana K t E Plg(M/M Q ) 



7ttV(7tt k + 7 t b E uIse ) and 



7t b E ulg 7(7t d E Sk + 7t b E Ulge )- Th e fractional contributions K tE {t { ° } ) for 
the chosen values of are also listed in TableQ If the uncertainty 
in tE is less than 20 per cent, it does not have a significant effect on 
the probability density. 

For the previously chosen selected values of t^ , 
(lg(M / Mq)^ and <t\ s {m/Mq) and as well as their expo- 
nentiated values are listed in Table |2| while the top panel of 
Fig. |H shows these values as a function of tE- While a mass 
M ~ 0.36 Mq for t^ ' = 20 d is in rough agreement with 
estimates using a 'typical' fractional lens distance x and transverse 
velocity v, the assumed mass spectrum with a low abundance 



4.3 Effective velocity and Einstein radius 

By eliminating x by means of the delta-function, Eq. IA15i yields 
for the probability density of the velocity parameter ( for a fixed 



an 



P?\C,vl E ) ) = 



r„ 



47, fE (% ( E ) ) 



C 4 x 



t'u 



6(mo - C 2 ) 



Mq*"" 

$>e (c,i(l±Vl-C7Mo 
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^Q^Vwv^)) d M0 



(0) 



(17) 
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Table 2. Expectation value and standard deviation for \g(M / Mq), x, and lg ( for selected values of the event time-scale t^. 



t B /(ld) 


(lg(M/M )) 


°lg (M/M ) 


M/Mq 




(x) 




<igC> 


°"lgC 


£/(100 kms- 1 ) 


/ V 


5 


-0.95 


0.39 


0.11 


2.4 


0.88 


0.10 


0.42 


0.15 


260 


1.4 


10 


-0.70 


0.37 


0.20 


2.3 


0.83 


0.12 


0.33 


0.16 


216 


1.4 


20 


-0.45 


0.39 


0.36 


2.5 


0.77 


0.15 


0.21 


0.17 


161 


1.5 


40 


-0.18 


0.46 


0.65 


2.9 


0.69 


0.20 


0.07 


0.21 


117 


1.6 


80 


-0.12 


0.58 


1.31 


3.8 


0.63 


0.23 


-0.07 


0.26 


84 


1.8 



In addition to the expectation value (lg(M / Mq)j and the standard deviation o"i g (M/M \, the corresponding exponentiated values 

M/M & = exp lo [(lg(M/M )}] and S M/J ft = exp 10 [<T lg(M / MQ )] are listed. Similarly, v = exp 10 [(lgC>] v c and = exp 10 [a lsC ] are given. The 
source has been assumed to reside in the Galactic bulge at a distance Dg = 8.5 kpc in the direction of Baade's window (I, b) = (1°, —3.9°) and the Galaxy 
model described in AppendixlBlhas been adopted. None of the listed values changes significantly if a 20 per cent uncertainty in is considered, where the 
distributions widen by less than 2 per cent, while mass and velocity estimate shift by less than 0.7 per cent, and the fractional distance x shifts by less than 

0.002. 
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Figure 4. Probability density of the fractional lens distance x for selected 
values of the observed event time-scale t&. Similar to Fig. |3}, the individual 
contributions Kt E p% of disk and bulge lenses are shown together with 
the total probability density p x , where the relative weight of the two lens 
populations is listed in Tablefflfor the chosen values of t e . 
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Figure 5. Expectation value and standard deviation for the logarithmic lens 
mass lg(M/Mpj), the fractional lens distance x and the logarithmic veloc- 
ity lg f as function of the event time-scale t e . While the solid lines mark 
the expectation values, the dashed lines limit intervals corresponding to the 
standard deviation. 
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Figure 6. Probability density of the velocity for selected values of the ob- 
served event time-scale £e- As for Figs. |3] and |4] the individual contribu- 
tions Kt E Pigf of disk and bulge lenses, with the corresponding weight 
factors listed in TablelTI are shown together with the total probability den- 
sity pig £ of lg where £ = i)/v c is the dimensionless velocity parameter 
and v c = 100 kms - 1 has been adopted. 



With X = yjl - ClJiQ, one also finds equivalently 



r„ 



4> a 



c 



(0). 



C 2 e( M r x -C 2 ) x 
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(0) 



x Z)*c(c,i(i±i)) *.(i(i±*)) di, 
± 

where the integration limits are given by 



(18) 



y/i - CVaC" for C<Po nin 
for C^Mo nin ' 



VWVmT for C<Mo 
for C > Mo 



max 


max 



(19) 



follows that of the velocity v for any value of the event time-scale 
£e- More precisely, if one defines te,o = Vt^ r E as the Einstein 
radius of the 'typical' mass Mo, corresponding to v = v c and x = 
0.5, one finds that pe = r^/r^fi = (, so fhatp PE (pE) = f>c(0- 

As for the lens mass M and the fractional lens distance x = 
Dh/Ds, expectation values and standard deviations for the trans- 
verse velocity v — Dl p = C v c at the lens distance are shown 
in Table [2] for selected time-scales is, whereas Fig. |6| shows the 
corresponding probability densities. As before, the source has been 
assumed to be located in the direction of Baade's window, at a dis- 
tance Ds = 8.5 kpc. The expectation value of lg£ as well as its 
uncertainty as function of the time-scale t& is also displayed in the 
lower panel of Fig. [5] 



5 CONSTRAINTS FROM PARALLAX AND 
FINITE-SOURCE EFFECTS 

5.1 Model parameters providing further information 

Further information about the lens mass M, the fractional lens dis- 
tance x, and the effective velocity parameter £ exceeding that pro- 
vided by the event time-scale £e can be obtained from events whose 
light curves are significantly affected by either the annual Earth's 
motion around the Sun or the finite size of the observed source star 
or even both of these effects. Either of these provides an additional 
relation between (M, x, £) from a model parameter that relates the 
Einstein radius te{M, x,Q to a physical scale which is either the 
Earth's orbital radius of 1 AU or the radius ii* of the source star. 

Moreover, for a binary lens, the total mass arises from the pe- 
riod P a nd the semi-major axis a, so that, as discussed by D ominikl 
( 1998b), further information about the lens properties arises from 
the lens orbital motion. However, it is quite difficult to obtain re- 
liable measurements of the full set of orbital elements in order to 
determine the period P and the parameter p = a/rs, which would 
provide a relation between M, x, and f . As pointe d out in the 
discuss ion of the event MACHO 1997-BLG-41 by Albr ow et all 
( 2000a), the lowest-order effects can be attributed to the actual pro- 
jected differential velocity between the components, which restricts 
only a subspace with two measured model parameters, while leav- 
ing another three undetermined. This strongly limits the power to 
constrain the lens and source properties. A proper discussion would 
be quite sophisticated and needs to be tailored to specific cases, so 
that it exceeds the scope of this paper. 

If the light curve is significantly affected by annual parallax 
resulting from the revolution of the Earth around the Sun, one can 
determine 7te = tvls /fe as a model parameter. In analogy to the 
definition of rj^ by Eq. <1 II . a corresponding dimensionless pa- 
rameter reads 



2tt 



(0), 



r E,0 



1 AU 



2.85 7T 



(0) 



1 kpc 



1/2 



(20) 



Since te = £e v, the distribution of the Einstein radius r& 



Similarly, if the finite size of the source star has a significant 
effect on the microlensing light curve, the time-scale t* = 
£b [R+/(Ds fe)], in which the source moves by its own angular 
radius relative to the lens, can be determined as additional model 
parameter. Alternatively, one might use a source size parameter 
p* — £*/£e instead. A corresponding dimensionless parameter can 
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Table 3. Fractional contributions K? lsl i and ft? ul f c of disk or bulge 
lenses to the event rate density yt E j7I - E (*e , tte ) for a Galactic bulge source 
by lenses in the Galactic disk or bulge, respectively, for a typical and 
selected 7Te ■ 







^.disk 


^bulge 

iE^E 


20 




0.35 


0.65 




0.015 


0.05 


0.95 




0.06 


0.19 


0.81 




0.25 


0.80 


0.20 




1 


1.0 


3 x 10~ 6 



With the event rate density 7t E ,7r E given by Eq. K^isk 



be defined as 



2tl O) r E , 



613 p^Q 



Ds 
1 kpc 



1/2 



(21) 



with the convenient abbreviation p^ q = {t*/t^) (R*/Rq) , 
where i?Q w 6.96 x 10 5 km denotes the solar radius. 
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Figure 7. Bivariate distribution of the event rate with the time-scale £e and 
the parallax parameter 7Te = ttls /£>e ■ For a bulge source at Ds =8.5 kpc 
in the direction of Baade's window (I, b) = (1° , —3.9°) and the Galaxy 
model described in Appendix|B] contours of 7ig[t E /(i[j)],lg7r E are shown 
at the levels 0.025, 0.05, 0.1, 0.2, 0.4, 0.8, and 1.6. 



5.2 Parallax 

From the expression for the event rate density for an observed 
and the related rjl E as given by Eq. 1121 . one finds with the addi- 
tional constraint S[ri^2 - a/(1 - x)/x (M/Mq)~ 1/2 ] the event 
rate density in both r)fj and rf^ E to be 

TV B ,17 



*E ' ^E / 
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X $ 



(M/M ) 3 



[i + [ v ^] 2 (m/m q )Y 
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x $> x 



1 + [fjffi] 3 (M/Mq) ' 1 + (M/M ) > 

1 



1 + (M/M©) 



d(M/M ) : 



(22) 



while 

/i(0) (0)\ 

2tj c 



1 AU 



7*7* 



E,0 



•i (0) 2 



T'E, (0) 
1AU %E 



(23) 



Hence, the bivariate distribution of the time-scale t E and the paral- 
lax parameter 7te is given by 7t E ,ir E = 7t E ,ir E /r, which is shown 
in Fig. [7] 
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Figure 8. Fractional contribution re? lsl i of disk lenses and k?" 1 ? 6 of 



bulge lenses to the total event rate density 7t E 



as a function of the 



time-scale and the parallax parameter 7Te for a Galactic bulge source 
at Dg = 8.5 kpc towards Baade's window. The set of contour levels at 
0.001, 0.01, 0.05, 0.1, 0.2, 0.5, 0.8, 0.9, 0.95, 0.99, and 0.999 corresponds 

,-disk 1 ..bulge 



to both quantities, since k? 1sj 1 +n'f _ 



1. The contour at 0.5 is shown 
in bold, while the contours at 0.1 and 0.9~are shown in light bold. 



values. The contributions of the Galactic disk given by «:t E * E 



-ydisk /(-disk , bulge 
'^e.^e/ \ '^E^E 1 n E ,ir E 
bulge 



tion k 



*E ! 7r E 



= It 



bulge 



), while bulge lenses contribute the frac- 
/(7 t d E *E + 7&£), are shown in Fi g-IH1 

(0) 



The measurement of the parallax parameter ir E in addition The corresponding values for the typical time-scale t E = 20 d 



to the event time-scale alters the fractional contributions of 
the individual lens populations to the event rate density for these 



and a few different 7r E are also listed in Table [3] Again, a bulge 
source at Ds = 8.5 kpc in the direction of Baade's window has 
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been assumed. By considering the contours at 0. 1 and 0.9 together 
with the distribution of tE and 7te as shown in Fig.0 one sees that 
strong preferences for one of the populations are not unlikely to be 
provided, whereas Fig. |2| shows that from £e the maximal prefer- 
ence for any time-scale is ~ 0.8 in favour of the bulge, achieved for 
ts, ~ 8 d. For tE > 2.5 d, smaller values of tte favour the lens to 
reside in the Galactic bulge, while larger tve favour the disk as lens 
population. For smaller time-scales, there is an intermediate region 
where this order is reversed. 

From the expression for the event rate density 7^ tV „ E for an 

event with given 4 0) and n^ 1 , given by Eq. <22> . one immediately 
finds the corresponding probability density of fio = M /Mo as 

8r r__(0)!5 ^(0)14 
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fn (0) v (0) ) 

; Vlt E ' '/»E I 
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X 4>f 
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o (°) (°) 
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(0) (0) 



1 + Mo 



(0) (0) 



1 + Mo 



(o) (0) 



(24) 



The probability density of the fractional lens distance x for ob- 
served tg 3 ' and 7Tg ' can be obtained with Eq. <1 5 1 by applying the 
parallax constraint in the form S[rj^ — yj (1 — x) /(mo X )/Vt°^]> 
yielding 

8r 1 



7vt E ,v* E 

x (l-x) 3 $ M/MQ ([<>] 
2(1 -x) 



(0) vi°2) 
-2 1- 



(0), 



X $ 



C 1 « (0 V 0) 

'/t E 'Ke 



#«(ar) 



For deriving p^^C; %g , »7we), one can first eliminate a; with the 
time-scale constraint in Eq. JA15I in order to obtain Eq. 1171 
and then use the parallax constraint in the form 5[rj^2 - (1 T 
\Jl — C 2 /mo)/( j 74e' ) C)]> or alternatively first eliminate /xo and 
then use 5[r))^ — 2 (1 — x )/(Cvi^)] as parallax constraint. Re- 
gardless of the way of approach, one obtains 

„(<>)//. „(o) „(o)s _ Tj 1 „(") 



(25) 



ft* 



x c e 1 



A^tO) (0) 



X $ 



(0) 



M/M 15) 



(26) 



For a typical time-scale = 20 d and different values for 
the parallax parameter tt^\ the probability densities P\ & (m/Mq), 
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Figure 9. Probability densities of lg(Af/MQ ), the fractional lens distance 
x and lg£, where £ = v/v c and v c = 100 kms -1 , for tE = 20 d 
and tte = 0.015, 0.06, 0.25, or 1.0 (solid lines) as well as that for an 
uncertain tte, i.e. based solely on the time-scale tE (dotted line). The source 
is thereby located in the Galactic bulge at Dg = 8.5 kpc in the direction 
of Baade's window (I, b) = (1°, -3.9°). 



p x , and pi g( ; are shown in Fig. [9] whereas expectation values 
and standard deviations for the related quantities are listed in Ta- 
ble |4] For comparison, the previously obtained results for an un- 
certain 7Te, i.e. based solely on the measured time-scale t£' are 
also shown. In most cases, the measurement of the parallax pa- 
rameter results in a significant reduction of the width of the distri- 
bution, equivalent to a reduction of the uncertainty of the consid- 
ered lens property, where the mass estimate improves most signifi- 
cantly. If, however, the parallax constraint forces the lens properties 
to fall into an a-priori disfavoured region, the expectation value is 
strongly shifted and the distribution may widen. The uncertainty is 
still dominated by the mass spectrum, mass density and velocity 
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Table 4. Expectation value and standard deviation for the lens properties \g(M / Mq), x 


and l£ 


£ for an 


event with time-scale *e and parallax 


parameter 7Te 
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0.20 


219 1.6 
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1.5 
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0.05 


0.28 


0.16 


190 1.4 




0.25 


-0.85 
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0.12 


129 1.3 




1 


-1.29 


0.34 


0.05 


2.2 


0.25 


0.01 


-0.19 


0.08 


64 1.2 



The row with 7Te marked ' — ' corresponds to an uncertain parallax parameter, i.e. the estimate is based solely on trj. Also listed are the exponentiated values 
M/Mq = exp lo [(lg(M/M ))] and? M/A -y = exp 10 [o- lg(M/M } ] as well as v = expi [(lgC>] «c and?„ /c = exp 10 [cr lgC ]. The source is located in 
the Galactic bulge at Dg =8.5 kpc in the direction of Baade's window (l,b) = (1° , — 3.9° ), and the Galaxy model described in AppendixlBlhas been 

adopted. 



distributions as compared to the contribution arising from the finite 
width of the time-scale distribution. As a result of the finite lim- 
its on the lens mass from the spectrum $m/m and the condition 
^ x ^ 1, the probability densities of the lens properties may face 
sudden cut-offs. 



5.3 Finite source size 

For a finite-source event with observed t 



and t+ and the re- 



lated rjl^ and rjj^ as given by Eqs. jilt and <21t . the event rate 
density Jr,^,^ (r^ >»7i° ) results from Eq. {T2j by applying the 
additional constraint 5\r]fJ - y/x/(l - x) (M/M )~ 1/2 ]. If one 
compares this with the case of parallax effects, one finds that r)^2 
assumes the role of tjI while x and 1 — x are interchanged, which 
is reflected in the result 
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The event rate density in (<e, t*) follows directly as 
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so that with p t = (t^/tv) (R*/Rq) , one finds (£e, p* ) to 
follow the distribution 
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Figure 10. Bivariate distribution of the event rate with the time-scale 
*E and the finite-source parameter q = (i,*/trj) (-R*/-Rq) _1 - The 
source is located in the Galactic bulge at Dg = 8.5 kpc in the direc- 
tion of Baade's window (I, b) = (1°, —3.9°). The plot shows contours of 
7lg[t E /(id)].l g 7r E correspond to the levels 0.025, 0.05, 0.1, 0.2, 0.4, 0.8, 
and 1.6. 



Fie. II II shows the fractional contribution ftt E ,t„ to the event rate 
density for a given £e and t* for the lens residing in either the 
Galactic bulge or disk. The measurement of finite-source effects 
turns out to be less powerful than that of parallax with most of 
the likely values not providing strong preference for either of the 
lens populations. However, small t+/tE Bulge lenses are preferred 
for intermediate values 0.0015 < {U /t E )(R* / Rq)' 1 < 0.007, 
where a strong preference however can only arise for tE < 10 d. 
Measurements of a small t*/t-E can provide a very strong prefer- 
ence for the lens to reside in the disk. 

From Eq. 1271 . the probability density of po follows as 
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Figure 11. Fractional contribution nf^ ^ of disk lenses and «^ lge of 
bulge lenses to the total event rate density 7t E ,t„ as a function of the time- 
scale and the finite-source time-scale t* = t E [R*/(Dg 6 E )] for a 
Galactic bulge source at Dg = 8.5 kpc towards Baade's window. With 
nf iai L + re° ulge = 1, the set of contour levels at 0.001, 0.01, 0.05, 0.1, 



0.2, 0.5, 0.8, 0.9, 0.95, 0.99, and 0.999 corresponds to both quantities. The 
contour at the level 0.5 is shown in bold, while light bold has been used for 
the contours at 0. 1 and 0.9. 



Table 5. Fractional contributions k^ 1s1 J and re bulgc of disk or bulge lenses 
to the event rate density 7t E ,t* (i El **) for a Galactic bulge source by 
lenses in the Galactic disk or bulge, respectively, for a typical t E and se- 
lected p i0 = (t*/t E )(iV%) _1 - 
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(30) 



while for the probability density of x, one finds 
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After elimination of /in using the constraint provided by 

the finite-source constraint becomes S[rj^ — 2 a;/(C )], while 

the elimination of x yields the constraint 6[rjl°' — (1 ± 

\/l - CVA*o)/('7tE C)]> so that either with Eq. (T7J or directly 
from Eq. IA15i . one obtains the probability density of £ for mea- 
sured tg ' and ti°' as 
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(32) 



Fig. 1121 shows the probability densities of \g(M / Mq), the 
fractional lens distance x, and lg £ for a microlensing event on a 
bulge source at Ds ~ 8.5 kpc in the direction of Baade's window 
for which t E and t* have been determined, where a typical time- 
scale = 20 d has been assumed, whereas a few different val- 
ues for p'q have been chosen. For the same values, Table|5|shows 
the fractional contributions of the bulge and disk lenses to the event 
rate density 7t E ,t* (t^ , t* " 1 ). While the parallax measurement has 
been found to provide the most effective reduction of uncertainty 
for the lens mass, one finds that the finite-source parameter p t q or 
the related time-scale £* most significantly affects the uncertainty 
in the effective transverse velocity v — v c £ or the related Einstein 
radius r E = t-E v. Some distributions show two peaks correspond- 
ing to the bulge and disk population. 



5.4 Combination of parallax and finite-source effects 

If both 7te and t* are measured, the lens mass M, its fractional 
distance x, and the effective velocity v are determined, so that 
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(33) 



and the uncertainty in these quantities is solely given by the dis- 
tributions of the model parameters 4b, tte, and i* or the related 
dimensionless rjt E , r] 7TE , and r\t t , respectively. 

As result of a fundamental property of logarithms and the lin- 
earity of the expectation value, the expectation value of the loga- 
rithm of a product of arbitrary powers of quantities £j separates into 
the sum of multiples of the expectation values of the logarithms of 
the individual quantities, i.e. 



(34) 



Similarly, one finds for the variances 
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Table 6. Expectation value and standard deviation for the lens properties lg(M/M ), x, and lg £ for an event with time-scale and finite-source 
parameter p^ Q = (**/*e) (R^/Rq)^ 1 . 
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M/Mq 




(x> 


fx 


<igC> 


°"lgC 


■0/(100 kms" 1 ) 


V / V 


20 




-0.45 


0.39 


0.36 


2.5 


0.77 


0.15 


0.21 


0.17 


161 


1.5 




0.0005 


-0.71 


0.46 


0.19 


2.9 


0.18 


0.18 


0.00 


0.33 


99 


2.2 




0.00125 


-0.26 


0.36 


0.56 


2.3 


0.71 


0.16 


0.35 


0.12 


222 


1.3 




0.003 


-0.67 


0.33 


0.22 


2.1 


0.84 


0.11 


0.04 


0.07 


111 


1.2 




0.0075 


-1.05 


0.43 


0.09 


2.7 


0.92 


0.08 


-0.31 


0.04 


49 


1.1 



For the row with the entry ' — ' for p t q, the estimate is based solely on 4e, while the finite-source parameter has been considered as uncertain. In addition to 
the basic estimates, the exponentiated values M/Mq = exp 10 [(lg(M / Mq)^] and ^j^/M = ex Pi0 [ CT lg(Af/AfQ)] as we " as ti = exp 10 [(lg £)] v c and 
S v /$ = exp 10 [(Ti g ^] are listed. The source is located in the Galactic bulge at Dg = 8.5 kpc in the direction of Baade's window (I, b) = (1°, —3.9°), and 

the Galaxy model described in Appendix iBlhas been adopted. 
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where Cov(xi, Xj) — Cov(xj,Xi) denotes the covariance of the 
quantities Xi and Xj, and Cov(xi, xi) = Vax(xi). 
While one naively finds the lens mass as 



M = 



c 1 1 AU t E 



7TE U 

taking into account the finite uncertainties yields 
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and with 8+ = R+/Ds not being correlated with the model param- 
eters tE, £*, and 7te, one obtains 
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+ Var (lg tt e ) + Var (lg + Cov (lg ^ 
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5.5 Limits arising from the absence of anomalous effects 

Frequently, anomalous effects such as those caused by the an- 
nual parallax or the finite source size escape detection from the 
photometric data taken over the course of the microlensing event. 
However, the absence of significant deviations from a lightcurve 
that is compatible with an ordinary event places upper limits on 
the model parameters 7te or t*. Rather than "defining" a cer- 
tain value by means of (^-functions, these constraints can be in- 
corporated by including O-functions in the respective expressions 
for the probability and event rate densities. With r/ 7rE = [(1 — 
x)/x] 1/2 (M/M Q )- 1/2 and jfr B < ry™ E ax , one finds a lower limit 
on the fractional lens distance 

x > x min = \ (39) 

i + (^r) (m/mq) 

for a given mass M. In analogy, for the finite source size, one finds 
with 774* = [x/ (1 
upper limit 



X < X 



)] 1 / 2 {M/Mq)- 1 ' 2 and r, u < vZ 
{r^Trf {M/Mq) 



C (M/Mq) 



(40) 



Table 7. Constraint on the fractional lens distance x = D^/Dg arising 
from upper limits on the annual parallax or the source size. 
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Vtt E 




^min 




0.1 Mq 


0.2 Mq 


0.4 Mq 


0.8 Mq 


0.015 


0.12 


0.998 


0.997 


0.994 


0.988 


0.06 


0.50 


0.98 


0.95 


0.91 


0.83 


0.25 


2.1 


0.70 


0.54 


0.37 


0.22 


1 


8.3 


0.13 


0.07 


0.03 


0.02 






vn 




X 


max 




P*,Q 


0.1 Mq 


0.2 Mq 


0.4 Mq 


0.8 Mq 


0.0005 


0.89 


0.07 


0.14 


0.24 


0.39 


0.00125 


2.2 


0.33 


0.50 


0.67 


0.80 


0.003 


5.4 


0.74 


0.85 


0.92 


0.96 


0.0075 


13 


0.95 


0.97 


0.986 


0.993 



The source star has been assumed to be located at a distance D$ = 8.5 kpc 
in the direction of Baade's window, = 2[r E q/(1 AU)]7te, 7?t* = 
2 (»E,©/%)P*,©- 



Taking into account these limits yields e.g. the probability density 
of the lens mass fio = M/Mq in generalization of Eq. I13i as 



pfoiWH^) = T775T7 V^o *Af/A/ M I Mo 



4 To 



x (mo) 



$ C (2 y/flo X(l - xj,Xj X 

S m,n (Mo) 

x [x(l-x)] s/3 & x (x)dx. 



(41) 



For a few selected masses, the resulting constraint on the frac- 
tional lens distance x = Z?l / Ds that arises for selected limits for 
the parallax or the source size is shown in Table0 where the 'stan- 
dard' source at Ds = 8.5 kpc in the direction of Baade's window 
has been assumed. Both the parallax and the finite-source constraint 
more strongly restrict smaller lens masses, while larger masses re- 
main possible at small distances with the parallax limit and at large 
distances with the finite-source limit. For the same parallax and 
source-size limits as listed in Table0 Fig. H3l shows the resulting 
probability density of the lens mass assuming an event with time- 
scale tE = 20 d for a source located at Ds = 8.5 kpc in the 
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Figure 12. Probability densities oi\g(M / Mq), the fractional lens distance 
x and lg£, where CJ = v/v c and v c = 100 kms -1 , for = 20 d and 
p^ = (U/IeX^/Bq)" 1 = 0.0005, 0.00125, 0.003, or 0.0075 
(solid lines) as well as that for an uncertain t*, i.e. based solely on the 
time-scale Je (dotted line). The source is located in the Galactic bulge at 
.Dg = 8.5 kpc in the direction of Baade's window (I, b) = (1°, —3.9°). 



direction of Baade's window (I, b) = (1°, —3.9°). Significant dif- 
ferences arise for tte < 0.8 or p^ q < 0.005. 

For the annual parallax, the transition between a geocentric 
and a heliocentric coordinate system does not influence the light 
curve and the orbital velocity is effectively absorbed into the event 
time-scale tE by contributing to the effective absolute perpendic- 
ular velocity. Therefore, it is the acceleration of the Earth 's orbi t 
that produces the lowest-order deviation (e.g. Smith et al. 2003). 
Within tE, this acceleration induces an angular positional shift of 
27T 2 7tls [*e/(1 yr)] 2 . so that k* = 27t 2 7r E [te/(l yr)] 2 is a suit- 
able indicator for the prominence of parallax effects. For an event 
time-scale £e = 20 d, a limit tte ^ 0.8 can be detected with a sen- 



D 5 



8.5 kpc, (i,b) = (l°,-3.9°) 




M/M 

Figure 13. Probability densities of \g(M / Mq ) for an event with time-scale 
tE = 20 d for a source located at Dg = 8.5 kpc in the direction of Baade's 
window (l,b) = (1°,— 3.9°), where an upper limit for the parallax 7te 
or the source size parameter q = (t*/tE) (R* / Rq)~ 1 applies. The 
curves shown as dashed lines correspond to the limit-free case. In the lower 
panel, this curve practically coincides with that for p 4 q = 0.0075. 



sitivity to ~ 0.05, whereas k w = 1 is reached for tE ~ 90 d, so 
that much smaller parallax limits can be obtained from such long 
events. If lens binarity can be neglected, finite-source effects be- 
come apparent if the angular source size 6 t becomes a fair fraction 
of the angular impact uo #e between lens and source. By requir- 
ing u < 2(0*/fe) = 2p„ q (R-^/Rq), a limit p+ < 0.005 
for Ri, — Rq is detected if uq < 0.01, corresponding to a peak 
magnification Ao > 100, whereas an impact parameter uq < 0.1, 
corresponding to Aq > 10, is sufficient for R* = 10 Rq . 



6 ESTIMATING ANOMALY MODEL PARAMETERS 

In order to judge whether any anomaly is likely to have a significant 
effect on the light curve, it is useful to estimate the value of param- 
eters that quantify the considered anomaly. As already pointed out 
in Sect. [5] the size of parallax effects arising from the orbital mo- 
tion of the Earth can be modelled by the parameter 7te = vtls/^e- 
With 



1 AU 



T E,0 



2r 



e,0 



(42) 



being the value that corresponds to a solar-mass lens at x = 0.5, 



one can define a 'typical' parallax parameter 7Te,o = [Vt^] 1 ^e 
for a given t^ ' and the chosen v c , where rj^ is defined by Eq. <1U . 

The corresponding ratio f}K E = 7te /tte,o is related to the basic 
properties as r)^ E = [(1 — a;) /x] p 1//2 , so that Eq. IA15I yields 
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Figure 14. Expectation values (solid lines) of lg 7te (top panel) and 
IgP* = 1§[(**Ae) (R* / 1 ] as a function of the event time- 
scale 4e for a source in the direction of Baade's window located at Ds = 
8.5 kpc. The dashed lines indicate limits defined by the standard deviation. 



the probability density of fj nE as 



(o) , - (o)\ _ 8rp „ 4 f Mo 

2 N 



(0) 



x 4>r- 



x <b x 



2 Mo 1 
1 + Mo i?i E ' 1 + Mo ^ E 

dp • 



Mo Vtte 



(43) 



From the respective definitions, one finds that Pfj^ E (v^e I *7t E ) = 

fae ) ^e/^e^/^e^e We)]- with 7tk b 8 iven b y 
Eq. fl2| and 7,^ ^ given by Eq. {22j. 

The top panel of Figure 1141 shows the expectation value of 
lg 7te along with its standard deviation as a function of the event 
time-scale £e- Since the variations in the basic system properties 
counterbalance each other with respect to the parallax parameter 
7Te, its expectation value shows only a slight variation with the 
event time-scale £e, while its variance is quite substantial. With 
K n — 27t 2 7TE [te/(l yr)] 2 being the angular positional shift in 
units of the angular Einstein radius 6 E induced by the acceleration 
of the Earth's orbit, which is a suitable indicator for the promi- 
nence of parallax effects, and 7te ~ 0.1 only weakly depending 
on the event time-scale, one approximately finds « w oc t B , where 
k,~6x for t E = 20 d, while k v ~ 0.1 for t E ~ 80 d. 

Finite-source effects can be studied in analogy to the par- 
allax case. With the definition fjt t = ^t, % E > so mat = 
[x/(l — x)] 1 / 2 Mo ly/2 > Eq. <A15> yields the corresponding prob- 



ability density as 



''0 



sro ^4 



Mo 



x $ 

X 



m/a/q [ Mo 
2 Mo Vt, 



(0) 



Mo Vt< 



1 + mo C 1 + A*o % 



X $ x | ; M0T?t *. o | d/i , 



1 + Mo f)t. 



where p,^ (j)*, ; 77™) = 7^,^, (% C e> ^/Vt^/ml^ntE (if)]' 
with 7 W given by Eq. 1121 and 7^ jr?ti given by Eq. J27t . As the 
bottom panel of Fieure fT4l shows, events with smaller time-scales 
are more likely to show prominent finite-source effects, where for 
a (giant) star with R ~ 10 Rq, one finds p* = i*/tE ~ 0.02 for 
t E = 20 d, whereas p* ~ 0.007 for i E = 80 d. 

If the lens that caused the microlensing event is a binary ob- 
ject and its orbital motion is neglected, its effects on the light curve 
are completely characterized by the mass ratio q, the separation pa- 
rameter d, where d# E is the angular instantaneous separation of 
its constituents (considered being constant during the duration of 
the event), and an angle describing the direction of the proper mo- 
tion with respect to their angular separation vector. The probability 
densities of the masses of the components Mi = M/(l + q) and 
M% = M [q/{l + q)] simply follow from the mass ratio q and 
the probability density of the total mass M as given by Eq. |l4j or 
by the corresponding relation given in Sect.[5]if parallax or finite- 
source effects are significant. An estimate for the instantaneous pro- 
jected physical separation r = d r E between the lens components 
can be obtained by multiplying the obtained d with the correspond- 
ing statistic for the Einstein radius te = D l 9e, so that 



(44) 



Pr 



(0) 



d(0) r E ,o 



,(°) 



(45) 



with pe = rE /rE,o, where te.o is given by Eq. 0, and pp E = 
defined by Eq. 1 17i . Beyond considering a fixed model parameter 
d, pf E can be convolved with its distribution. 

With x = f/o, denoting the projection factor between the 
semi-major axis a and the actual projected separation f , one finds 
a = r/x = (drri)/x. In addition to M, x, and the value of 
a therefore depends on the projection factor x as further prop- 
erty which is distributed as $x(x) as given by Eq. IC6L where 
^ X ^ X max - F° r the probability density, one therefore finds 



P™ E >/r E ,o;<V w ) 



(0) j(0K 



Pf/r EO {f/rE,o;v < t ° E \d {0) ) x 
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According to Kepler's third law, the orbital period reads P = 
2n [a 3 /(GM)] 1/2 , being a function of the total mass M and the 
semi-major axis a. With 



Po = 2n- 



'gg r E,Q 
GM 



(47) 



77P = -P/Po is related to the basic system properties as r)p = 

[2 (d/x)] 3 ^ 2 — x )] 3 ^ 4 Pc/ 4 > so m at the corresponding proba- 
bility density becomes 
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(49) 



Using the property for the expectation values and the variances 
stated by Eqs. 1341 and <35> . one finds that 

(lg(Mi/M )> = (lg(M/Mo)) - <lg(l + q)) , 

(lg(M 2 /M )> = (]g(Af/Afo)> + (lg q) - (lg(l + «)> , 

<lg(f/r E , )) = (lgp E ) + (lgd) , 

<lg(a/r E , )) = (lgp E ) + (lgd)-{lgX> , 

<lg(P/Po)) = ~ <lg(a/r E ,o)) - ~ <lg(M/M )> , 

= 5 [<lgPE) + (lgd)-(lg X )] - 

- 1 (lg(M/M )) , (50) 



and 

Var[lg(Mi/M )] 
Var[lg(M 2 /M )] 



Var [lg(M/M )] + Var [lg(l + q)] 
Var[lg(M/M )] + 

+ Var{lg[g/(1 + q)]} , 



Table 8. Estimates for physical properties of the two components of the 
lens system that caused microlensing event OGLE 2002-BLG-099 based on 
the model parameters reported by Jaroszvriski et al. 1 2004 ) and the Galaxy 
model described in AppendixlBl 



circular elliptic 


Mi [M ] 


0.38 (2.7) 


M 2 [M Q ] 


0.093 (2.7) 


D L [kpc] 


5.7 ± 1.7 


a [AU] 


6.2 (1.8) 5.8 (2.0) 


P[yr] 


22 (2.0) 20 (2.4) 



The estimates are based on the expectation values of x, \%{M/Mq), lg p E , 
and lgx, where a source distance of Dg = 8.5 kpc has been adopted. 
Numbers in brackets denote the uncertainty factor that corresponds to the 
standard deviation of the logarithm of the respective quantity. Afi and M 2 
are the masses of the primary and secondary component of the lens binary, 
respectively, Dl is their distance from the observer, a denotes the orbital 
semi-major axis and P denotes the orbital period. The latter quantities have 
been calculated both assuming circular orbits or elliptic orbits with the ec- 
centricity being distributed as <E> e = (4/tt) yl — e 2 . 



Var [lg(f/r E ,o)] 
Var [lg(a/r E ,o)] 

Var [lg(P/Po)] 



Var(lgp E ) + Var(lgd) , 

Var (lg p E ) + Var (lg d) + Var (lg X ) , 

9 

- [Var(lgp E ) + Var(lgd) + 

+ Var (lg x)] + ~ Var [lg(M/M )] - 

- jCov[lgpE,lg(M/M )] . (51) 

As an example, let us consider the microlensing event OGLE 
2002-BLG-099 Ijaroszvfiski et all2004h . where the observed light 
curve suggests the lens star to be a stellar or brown-dwarf binary, 
while the absence of observed caustic passages leaves the possibil- 
ity that the source is a binary instead. For the binary-lens model, 
the mass ratio is q ~ 0.25, while d — 1 .963 and £e = 34.4 d. 
No signals of annual parallax or finite source size have been de- 
tected, whereas more than 60% of the observed light at baseline 
arises from a source other than the lensed background star. 

With the Galaxy models described in AppendixlBl one finds 
the expectation values and uncertainties of the two components of 
the binary lens system, the distance from the observer, the semi- 
major axis of the planetary orbit, and the orbital period as listed 
in Table [5] while the probability densities of these quantities are 
shown in Figs. 1151 and 1161 For the deprojection of the orbit, ei- 
ther circular orbits or elliptical orbits, where the eccentricity is dis- 
tributed as <f> e = (4/tt.) Vl — e 2 , have been considered. More 
details about the statistics of binary orbits can be found in Ap- 
pendix Icl For circular orbits, (lgx) = —0.133, corresponding to 
a factor X = 0.736, so that a ~ 1.36 f. The standard deviation 
of aigx = 0.183 is equivalent to a factor 1.35, yielding an interval 
a ~ 1.01 . . . 1.83 f. In contrast, one finds for elliptic orbits with the 
adopted distribution of eccentricities an expectation value (lg x) = 
-0.104, which yields a factor \ = 0.787, so that a ~ 1.27 r. 
In this case, the standard deviation is (Ti g -% = 0.234, which corre- 
sponds to a factor 1.71, spanning an interval a ~ 0.74 . . . 2.18 r. 
The differences between circular orbits and the elliptical orbits ac- 
cording to the adopted eccentricity distribution, which are seen in 
the respective probability density of the orbital period, reflect the 
distribution of the projection factor < l > x(x) as shown in Figure IcTl 
While for circular objects, a dominant contribution results from 
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Figure 15. Probability densities of the mass Mi of the lens star and its 
planet M2 as well as for their distance Dl from the observer for the 
microlensing event OGLE 2005-BLG-071 and the binary-lens model re- 
ported by Jaroszvhski et al. 120041) . for which q ~ 0.25, d = 1.963 and 
t E = 34.4 d. 



X < 1, the adopted elliptical orbits produce a small excess for large 
orbital periods and a significant tail towards smaller orbital periods 
due to projection factors 1 < \ < 2. 



7 DETECTION EFFICIENCY MAPS FOR LENS 
COMPANIONS 

As pointed out bv iMao & Paczvnskil Il99ll) . the light curve of a 
galactic microlensing event may reveal the existence of compan- 
ions to the lens star, which includes stellar binaries as well as plan- 
etary systems. If the orbital motion of a binary lens does not result 
in a significant effect, only two characteristics of the binary influ- 
ence the light curve, namely its mass ratio q and the separation 
parameter d = f/r-E, where r is the actual projected separation 
perpendicular to the line-of-sight. The detection efficiency e(d, q) 
for a companion to the lens star with a separation parameter d and 
a mass ratio q is defined as the probability that a detectable signal 
(event S) would arise if such a companion exists (event C(d, q)), 
i.e. 

e{d,q)=P(S\C(d,q)) . (52) 

For a given event, a de tection efficiency map (e.g. lAlbrow et alJ 
2000b : lRhie et alJl200d) can be calculated for any combination of 
the parameters (d,q). Let F(d~ , d + , q~ , q + ) denote the average 



OGLE 2002-BLG-099 
D s = 8.5 kpc, (i,b) = (5.45°,2.98°) 



1.5 
1 

0.5 





1 1 1 


1 1 1 1 1 1 1 1 
circular - 




elliptic _ 




/ V 

/' V 

/ \ 

/ ' \ 






/ ' \ 
// v 

/ / 











10 

a/(l AU) 



100 




0.5 



10 100 

P/(l yr) 



Figure 16. Probability densities of the semi-major axis a (top) and the 
orbital period P (bottom) for the microlensing event OGLE 2002-BLG- 
099 IJaroszvriski et all2004l) assuming the reported binary-lens model, for 
which q ~ 0.25, d = 1.963 and i E = 34.4 d. The dashed lines 
correspond to the assumption of circular orbits, while curves drawn as 
solid lines correspond to elliptic orbits with an eccentricity distribution 

$ e = (4/ti) VT^. 



number of companions with distance parameter d £ [d ,d + ] and 
mass ratio q £ [q~ , q + ], so that 



d+ q+ 

F(d-,d + ,q-,q+)= J J f d , q (d,q) dddq. 

d~ q~ 

From a sample of TV events, one then expects to detect 



(53) 



H(d-,d+,q-,q + ) 



U, q (d,q)e { Z(d,q)dddq (54) 



d q 

companions, and the probability to detect at least one signal reads 
fc.f. lAlbrow et alfeOOll) 

i(d~ ,d + ,q~ ,q + ) 

N 



n 



fd, q {d,q) £ ( jUd,q) dddq 



(55) 



However, rather than obtaining information about the awkward 
fd,q{d, q), one would like to investigate the abundance of compan- 
ions (such as planets) as function of the physical properties such as 
the companion mass M2, the semi-major axis a, and the orbital 
eccentricity e. While microlensing does not provide a means to 
study the dependence on the orbital eccentricity, for which a dis- 
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tribution needs to be assumed, the adoption of a Galaxy model al- 
lows to compare the microlensing results with assumed abundances 

fa,M 2 (a, M 2 ). 

For a given projected separation f and a companion mass Ma, 
the probability density of (d,q) follows that of (te,Mi), so that 
the detection efficiency in these physical lens characteristics reads 



£f,M 2 {f,M 2 ) = 



£d, q (d, q) p d ,q{d, q; f, M 2 , rj^ ) dd dg 



cd., 



M 2 /Mq 



,<ve 



E,0 



//.() 



u u 

x p PE , MD (pe, Mo; % ( e) dpE dp , (56) 

where pe = r^/rEfi with te,o and Mo as defined in Sect. [4] 
Rather than to the total mass M, the time-scale hereby refers to 
the mass of the primary Mi = M/(l + q) and the mass spectrum 
is adopted as $ Mi/Mq (Mi/M ) = 6 (Mi - M 2 ) $a//m [(1 + 
q)Mi/M(7)]. Whereas for Mi — M 2 , one needs to distinguish be- 
tween close binaries, where the best-fit single-lens time-scale refers 
to the total mass, and wide binaries, where it refers to one of the 
constituents, for the relevant M 2 < 10 _2 ' 5 Mq discussed here, 
g< 1 and Mi sa M is a fair approximation. Moreover, a single 
t-E for each of the events, corresponding to the value estimated for 
a single lens, rather than an optimized t& for each pair [d, q) can 
be used, since as shown previously, shifts in is by less than 20 per 
cent can be safely neglected relative to the width of the broad dis- 
tributions of lens mass, distance, and velocity, and the uncertainties 
of the Galaxy models. 

The distribution of (pe, po) follows from Eq. 1171 as 

To 



p(o) 



(pE,Mo;»7t E ) = 
Pe Q(po - Pe) 

E 

± 



47, tF (% ( ° ) ) 



- Pe/po 



M/Mq 



po 



(0) 



4> 



X $ : 



(i(-vT 



Pl/Po 



(57) 



A given semi-major axis a of the binary lens may result in 
different projected actual separations r depending on the spatial 
orientation of the orbit, the orbital phase and the orbital eccentricity. 
With <E>^ (x) denoting the probability density of x = r/a as derived 
in Appendix|C| the detection efficiency in (a, M 2 ) follows as 



e a ,M 2 (a,M 2 )= J e f , M2 (xa,M 2 )$*(x)dx. (58) 
o 

For circular orbits, the expression for $x(x) given by Eq. <C7> 
yields with a variable substitution in favour of w — yjl — \ 2 and 

X max = 1 



a ,M 2 {a,M 2 ) = I £f,M 2 (y/l-w 2 a, M 2 )dw. 



(59) 



Follo wing a pilot a nalysis of the event OGLE 1998-BUL- 
14 lAlbrowet alfe oOOb). for which the underlying technique has 
been developed, the PLANET (Probing Lensing Anomalies NET- 
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Figure 17. Bivariate probability density Pi gr - E /(i AU),lg m/Mq 
tion of the Einstein radius rF, and the lens mass M for the time-scale 
t E = 97.4 as determined bv lYoo et alj <2004 for OGLE 2003-BLG-423, 
and the corresponding location (/, b) = (0.50°, —5.18°), while a source 
distance Dg = 11 kpc has been assumed. The bold line marks the upper 
limit for te and the mass cut-off at M = 0.01 Mq, which arises from the 
adopted mass spectrum. 



work) collaboration has calculated detection efficiency maps in 
the parameters (d, q) for its data collected from 1995 to 1999 
in order to derive upper abundanc e limit s on planetary compan- 
ions to the lens stars JAlbrow et alJl20oll Gaudi et alJl2002l) . De- 
tection efficiency maps have also been derived by the MPS (Mi- 
crolensing Planet Search) and MOA (Microlensing Observations 
in Astrophysics) coll aborations for the event MACHO 1998-BLG- 
35 (Rhie et al. 2000), and several other groups for suitable events 
fTsapras et alfeOOllBond et all2002l.lYoo et alJ2004[ | Dong eta?] 
2005^. while lTsapras et alJJ2003h and lSnodgrass et alji2004 have 
determined planetary abundance limits from OGLE (Optical Gravi- 
tational Lens Experiment) data. The largest sensitivity to planets so 
far was ac hieved for the ev ent MOA 2003-BLG-32/OGLE 2003- 
BLG-219 <Abe et aljEbol that showed an extreme peak magni- 
fication in excess of 500. PLANET is in the progress of carrying 
out a new comprehensive analysis including the more recently ob- 
served events JCassan e t al. 2005), where, based on the results pre- 
sented in this section, a planetary abundance f a ,M 2 (a, M 2 ) rather 
than fd, q {d, q) is considered. 

For the parameters of the event OGLE 2003-BLG-423, where 
the source star is located towards (l,b) = (0.50°, —5.18°) 
and the best-fitting event time-scale assuming a point lens is 
tE = 97.4 d, Figure fTTl shows the bivariate probability density 
Pigr E /(i AU),igM/A/Q as function of the Einstein radius te and 
the lens mass M, where Ds = 11 kpc has been assumed. The bold 
line marks the upper limit for te, which corresponds to the lens 
being half-way between source and observer (x = 0.5), and the 
mass cut-off at M = 0.01 Mq, resulting from the adopted mass 
spectrum. As the Einstein radius te approaches its maximal value, 
Pi g r E /(i AU),i g M/A/ diverges. 

The detection efficiency as f unction of t he mo del parameters 
(d, q) that has been calculated bv lYoo et alJ J2004I) based on data 
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detection efficiency e for lens companions 
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Figure 18. Detection efficiency maps resulting from MicroFUN and OGLE data collect ed for the event OGLE 2003-BLG-423. The top left panel shows 
the detection efficiency £d,q{d,q) in the model parameters (d,q) as calculated bv lYoo et alj 12004). where the actual angular separation between the lens 
constituents is d fig and q denotes their mass ratio. The remaining panels show the detection efficiency as function of the physical lens properties derived using 
the Galaxy model described in Appendix[E] using the event time-scale = 97.4 d resulting from the photometric data, D§ = 11 kpc and the Galactic 
coordinates (£, 6) = (0.50° , —5.18°). All these plots refer to the secondary mass M2, but different separations, where the top right panel shows the detection 
efficiency for the actual projected distance perpendicular to the line-of-sight r = dr^,, the bottom left planet refers to the orbital radius r assuming circular 
orbits, and the semi-major axis a is used in the bottom right panel assuming elliptical orbits with the eccentricity distribution <T> e = (4/tt) yl — e 2 (see 
Appendix let. For each of the plots, contour levels are shown at e = 0.1, 0.2, 0.33, 0.5, 0.67, 0.8, 0.9, 0.95, 0.98. While all these contours shown up for the 
plot of £d,q{d, q), some contours that correspond to larger detection efficiencies fall outside the displayed region for the other plots. 



collected by MicroFUN (Microlensing Follow-Up Network) and 
OGLE for this event is shown in the top left panel of Fig. 1181 Data 
from the OGLE survey are ma de available on-line 3 as the events 
are progressing lUdalski 2003), which significantly eases the as- 
sessment of their parameters and thereby allows the optimization 

3 |http : / /www . astrouw . edu . pi / ~ogle/ogle3/ews/ews .html| 



of follow-up observations. Using the expressions presented in this 
section and adopting the Galaxy model described in Appendix IrH 
corresponding detection efficiency maps in physical quantities have 
been determined, where the single time-scale £e = 97.4 d has been 
used rather than the more exact best-fitting value for each (d, q) that 
refers to the primary mass, and a source distance of Ds = 11 kpc 
has been assumed. These maps are shown in the remaining panels 
of Fig. 1181 where the detection efficiency refers to the secondary 
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mass M2 and the projected actual separation f , the orbital radius r 
for circular orbits, or the semi-major axis a for elliptical orbits. The 
distribution of the eccentricity e for elliptical orbits has been cho- 
sen to be & e (e) = (4/7t)Vl — e 2 , which provides a rough model 
of the eccentricities for planetary orbits found by radial velocity 
searches (see Appendixlcl, 

All detection efficiency maps show the similar pattern of a 
maximum efficiency for a characteristic separation decreasing both 
towards smaller and larger separations and a decrease of detection 
efficiency towards smaller masses. As compared to the detection 
efficiency in the model parameters (d, q), the other panels show the 
detection efficiency being stretched over a broader range of param- 
eter space, so that peak detection efficiencies are reduced, while 
smaller values occupy wider regions. The main broadening occurs 
on the transition from (d,q) to (f, M%), so that the width of the 
distributions of the lens mass, distance, and velocity yield the dom- 
inant contribution rather than the orbital projection, which has a 
more moderate but highly significant effect. While the map in the 
projected actual separation f reflects the upper limit of the Ein- 
stein radius r&, this is smeared out by the distribution of the projec- 
tion factor when considering the semi-major axis instead. The av- 
erage orbital radius r for circular orbits exceeds the average semi- 
major axis for the considered elliptic orbits, and for elliptic orbits, 
the detection efficiency is also stretched towards smaller a, with 
1 < X — r/a < 2 being possible. A substantial detection effi- 
ciency for small planetary masses results from the large abundance 
of parent stars with small masses, whereas stars much heavier than 
the Sun are rare, and the fact that large detection efficiencies from 
larger mass ratios provide substantial contributions with the finite 
width of the mass distribution for a given event time-scale. 



8 SUMMARY AND OUTLOOK 

This paper provides a comprehensive theoretical framework for the 
estimation of lens and source properties on the basis of the re- 
lated model parameters that are estimated from the observational 
data. This formalism can be used to answer a large variety of ques- 
tions about the nature of individual microlensing events. With the 
adopted Galaxy model and a source star residing in the direction 
of Baade's window (l,b) = (1°, -3.9°) at D s = 8.5 kpc, 35 
per cent of all ongoing events (not identical with the monitored 
or detected ones) are caused by lenses in the Galactic disk and 65 
per cent by lenses in the Galactic bulge. While the bulge lenses 
clearly dominate the total event rate only for very small time- 
scales £e < 2 d, the disk lenses yield the larger contribution for 
time-scales 2 d < t E < 40 d and t E > 100 d, where the lat- 
ter however yield only a small contribution to the total rate. For 
2 d < t E < 40 d, bulge and disk lenses yield comparable con- 
tributions. The provision of probability densities of the underlying 
characteristics of the lens and source system such as the lens mass 
M, the distance Dl and the effective transverse absolute velocity 
v under the assumption of mass spectra, mass densities and veloc- 
ity distributions yields the largest amount of information that can 
be extracted from the observations, i.e. much more than by a finite 
number of moments. While a mass M ~ 0.36 Mq for tE = 20 d 
is in rough agreement with estimates using a 'typical' fractional 
lens distance x and transverse velocity v, the assumed mass spec- 
trum with a low abundance for M > 1 Mq forces the expected 
mass to be more narrowly distributed with t E rather than to follow 
the naive M oc t| law. In particular, there are only 1.5 decades 
0.09 Mq ^ M ^ 3 Mq for the expected mass if time-scales 



3 d ^ *e ^ 150 d are considered where the inclusion of one stan- 
dard deviation extends this range to 0.03 Mq ^ M ^ 15 Mq. 

Additional constraints such as those resulting from a measure- 
ment of the relative proper motion between lens and source from 
observed finite-source effects or the relative lens-source parallax 
as well as upper limits on these quantities resulting from the ab- 
sence of related effects can be incorporated. Explicitly one sees 
how uncertainties in M, Dl, and v are reduced, although the re- 
spective probability densities can also widen if the additional con- 
straint forces the lens to assume values that fall into regions dis- 
favoured by the given time-scale. For any set of observables, one 
also obtains a probability for the lens to reside in any of the po- 
tential lens populations. Unless there are sufficient constraints to 
yield a sharp value for the lens mass, distance, and velocity for a 
given set of model parameters, the uncertainties of the latter can 
be neglected against the broad distributions of the relevant charac- 
teristics of the lens populations and the Galaxy model uncertain- 
ties. With significant effects by annual parallax on the light curve 
starting at 7te < 0.8, such a limit can be detected in an event 
with time-scale tE = 20 d with a sensitivity to an angular po- 
sitional shift within tE of K n 9e ~ 0.05 #e, whereas «„- = 1 
is reached for tE ~ 90 d. Similarly, finite-source effects become 
apparent if the angular source size becomes a fair fraction of 
the angular impact uo 9e between lens and source. By requiring 
uo < 2(0*/0 E ) = 2/O^q (7?*/7?q), a limit p+ Q sC 0.005 for 
R+ — Rq is detected if uo J5 0.01, corresponding to a peak mag- 
nification Aq > 100, whereas an impact parameter uo < 0.1, cor- 
responding to Ao > 10, is sufficient for i?* = 10 Rq. 

In addition to the basic quantities, probability densities of the 
orbital semi-major axis and the orbital period for binary lenses, as 
well as of any quantity that depends on the basic characteristics, 
can be obtained. The bivariate probability density of the Einstein 
radius ?~e and the lens mass M together with statistics of binary 
orbits yields detection efficiency maps for planetary companions to 
the lens star as function of the planet mass A/2 and its orbital semi- 
major axis a rather than of the model parameters d and q, where 
d#E is the actual angular separation from its parent star and q is 
the planet-to-star mass ratio. The presented formalism has been ap- 
plied to some first examples and will be used for discussing the im- 
plications of many further events. This paper explicitly shows the 
distributions of event properties for the binary-lens model of mi - 
crolensing event OGLE 2002-BLG-099 ijaroszvnski et alJl2004h . 
namely of the masses of the lens components and their distance 
from the observer, as well as of the orbital semi-major axis and pe- 
riod. Moreover, it shows the detection efficiency map in (a, A/2) 
resulting from MicroFUN and OGLE data <Yoo et all2 004) for the 
event OGLE 2003-BLG-423. As a function of (o, A/ 2 ), the detec- 
tion efficiency stretches over a much broader range of parameter 
space than for the (d, g)-map. In particular, this results in a larger 
detection efficiency for low-mass planets than one would expect 
from typical values. 
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APPENDIX A: GENERAL PROBABILISTIC APPROACH 

Let us consider a system characterized by k properties a,i (i — 
1 . . . k) that are distributed statistically, where Q ai («i • • • a k ) dat 
gives the probability to find the property ai in the interval [ai, ai + 
daj which might depend on all ai . ..ak- Further consider any re- 
alization of these system properties yielding a specific contribution 
to observed events described by a weight function Q(di . . . a k ) 
which may be chosen appropriately to include selection effects 
caused by the experiment, so that the total event rate according to 
their statistical representation is given by 

= Ji k \Jn(a 1 ...a k )^J[^ at (a 1 ...a k )^ j JJ day j(Al) 

where the notation refers to a fc-dimensional integral. 

Hence, the probability density of the properties ai . . .a k 
among all observed events is proportional to 



Pa x ...a h {ai ■ ■ ■ a k \ $ ai . . . <& akl fi) 



Ofai . . . a fe 



(A2) 



so that an appropriately normalized probability density is given by 
p ai ...a k (ai . . . a k \ $ ai . . . $ afc , fi) 

= p ai ...a k {ai . . .a k ;$ ai - ■ - ®a k ,&) j 

/ /•••/ p ai ...a k {ai ... a k ; $ ai ...$ afc ,fi) Y[da t , (A3) 

which does not depend on any constant factors in fi. This means 
that p ai ...a k (ai . . . a k ; $ 01 • • • $a fc , fi) is obtained by weighting 
the intrinsic statistical distribution Y\. =i ( a ») by f2(ai . . . a k ) 
and normalizing the resulting product, so that Eq. <A3> corresponds 
to Bayes' theorem. 

A specific event involves a set of n observed parameters 
/; (ai . . .a k ), where / = 1 . . . n, which in general depend on the 
k basic underlying properties ai, but are not necessarily identical 
to these. With specific realizations / ; ' ' for an observed event, the 
event rate can be written as integral over these realizations 

r($ ai ...$ afe) n) 

r(B) ^7/x.../»(/i (0) -/i 0) ;*-,. x 



(0) 



(A4) 



with the event rate density 

7/ 1 .../„(/i C0) ---^ 0) ;*a < --.*a fc! n) 



y^yWi . . .a*) |n* (/j(ai . . . ak ) - 



X 



(A5) 



so that the corresponding probability density of the basic properties 
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ai . . . a k is given by 

P ( al. ak {ax... a k ,fW . . . /< 0) ; $ 01 . . . <£v fc , SI) 

= [7/x.../»(/f 0) --^ 0) ;*. 4 -*«*.n)r 1 x 

x n(ai...ajk)| JJ* ryj(ai...o k )-/, (0) ) } x 
^ 1=1 > 

^ i=l ' 

while the probability density of a single property a r reads 
pM(a r ;rt K..ti 0) -,$ ai ...<S> ah ,n) 

= hA...f n (fi 0) ---fi 0) ;^ i ..^ ak ,si))- 1 x 

x J&rV y 0(oi . ..o h ) | JJ* (/,(ai . . . a k ) - /,<*'>) 

i=l > ^ 3=1 J 



(A6) 



(A7) 



If the observables /i . . . fn for one or more events follow a 
distribution i &f 1 ...f n (/i •••/»), the probability density of the ba- 
sic properties ai . . . ajb arises from an integral over the probability 
density p' ' for the fixed values /} ' . . . , given by Eq. jA6j, as 

Pai...a fc (oi---a*;*o 1 ...$a fe ,n) 

= jw[p£l. ah (ai...a k ,h...f n -,$ ai ...$ ah ,Sl) x 

n 

x */!.../» (A • •-/») II d/i - (A8) 

;=i 

If the observables are statistically independent, their distribution 
factorizes as 9 h ... fn (fi .../„) = JlLi */i (/*)• while fixed 
values of f correspond to the distribution \& f t (/; ) = 5(fi — f, ), 
distributions around a central value / ; ' ' with a standard deviation 
a /j can be approximated by the Gaussian distribution 



^Gauss (/i) 



2na fl 



exp ■ 



(A9) 



In case the n observables / = (/i,. ••,/«) are correlated, 
.../ n (A •■•/«) can t> e modelled as a multivariate Gaussian dis- 
tribution 

(A •-•/») = ' 



(27t)™/ 2 |C|V2 

exp{-i (/-/^""C-M/-/^)} 



(A10) 



where C 1 is the inverse and C| is the determinant of the covari- 
ance matrix C, and = (/f 5 ,...,f n 0) ). 

The moments of any property o(ai ... a^) for fixed values of 
the observables /} ' . . . /„ follow from the expectation values 

(/(ai...a fc )> (0) = |w/'/(ai...a,) x 

x p<°>. ..„» (ai . . . a fc , /< 0) . . . /< 0) ; * 01 ■ ■ ■ $a fc , Q) x 
fc 

x JJda,, (All) 

i=i 



and for a distribution f 1 ..j n , one finds in analogy to Eq. <A8> 
(/(ai...a fc )>* = |.Wy (/(oi...a fc )> (0) x 

n 

x */x.../„(a---a) (ai2) 



where interchanging the integrations over dai . . . da k and over 
d/i . . . df n yields the equivalent expression 



(/(ai...O)* = J. ':. J, r ■<, ...„,>: x 

X Pai...o fc (Ol • • • Ofc, A (0> • ■ • fn } ; $<>1 • ■ • $a fe , SI) X 

x ]Jda s . (A13) 



In particular, the standard deviation is given by 



o 3 = yj (g 2 (ai . . . a*)) - (g(ai . . . a k )Y 



(A14) 



Beyond the moments, one finds the complete probability den- 
sity of a general property g(a\ . . . a k ) for fixed values of the ob- 
servables f[ . . . fn to be 

pf (g;fi°K..f n V;3> ai ...* ak ,Sl) 

= [7/ 1 .../n(/i (0) ..-/J? ) ;*« 4 -*-*,n)r 1 x 

x Ji k) .Jsi(a 1 ...a k )8(g-g(a 1 ...a k )) x 

f[ 5 (fi(ai ...a k )- // 0) ) 1 1 Y[ *a 4 (ft) da, | ,(A15) 

1=1 > i=l J 

while for a distribution \P / x . . . / n , one obtains 
Pg(g;<f> ai ...$a k ,Sl) 

^-/pf A .../»; #„,... *« fc ,n) x 



(A16) 



*/i.../„ (A ■•■A) H i/ 



It is important to distinguish carefully the different quan- 
tities that have been defined in this section. The system 
properties oi . . . a k are distributed statistically among the 
population according to 5> ai . . . 5> a)j . With Sl(ai...a k ) be- 
ing the weight of any realization to the number of pro- 
duced events, one expects a\ . . . a k being distributed as 
p ai ...a k (oi . . . a k ; <E> ai . . . <E> afc , SI) among all events. For a given 
event, with a set of observables /i(ai . . . a k ) . . . / n (oi . . . a k ) 
being realized as f^ . . . f^ , one infers a stochastical proba- 
bility density pi? (oi . . . a k , f[° ' . . . fL 0) ; $ 01 . . . $ afc , fi) of 
oi...Ofc or pf'(g,fi . ../i 0) ;$ai . ..*a fc ,n) of any spe- 
cific property g r (ai...afc), which does not need to be an ob- 
servable /; or a basic property a,j. Finally, one can consider 
the observables /i . . . f n to follow a stochastical distribution 
for a single event and/or a statistical distribution from sev- 
eral events, namely ^fx...f n yielding the probability densities 
P*...o s (ai . . . o fc ; $ ai • • • $a fc , SI) or p*(p, $ ai . . . $ afc , fi). 
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Table Bl. Coefficients for the mass laws for different parts of the mass spectrum (following lchabriej <2003l) 1 
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0.33 
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For M m i n M $C Af max , either a power-law mass function ^i s (m/Mq) [^&(M/Mq)] oc (M/Mq) a 
oraGaussian distribution 4>i g(A f/M ) [lg(A//M )] oc exp{-0.5 [lg(Af/M ) - lg(M c /Af )] 2 /(cr c ) 2 } is 
adopted. 
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Figure Bl. Adopted mass function $\ & m/Mq 
lens objects in the Galactic disk or bulge whose parameters are listed in 
TablelETI 



APPENDIX B: MODEL OF THE GALAXY 
Bl Mass spectrum 

While the mass spectrum $>m(M) as defined by Eq. gives 
the decomposition of the mass density into objects with mass in 
the range [M, M + dM], the decomposition of the number den- 
sity is given by the mass functio n &m(M) oc M -1 $m(M). As 
pointed out e.g. in the review o f lChabrieJ 12003). the mass function 
&1 S (M/Mq) can be fairly well approximated piecewise by differ- 
ent kinds of analytic functions. In particular, for selected ranges 
M^j n ^ M ^ Mmax. one ma y consider a power-law mass spec- 
trum in M, i.e. 



(i) r 
5(M/Mq)L 



i(M/M Q )] <x {M/MqY 



(Bl) 



with the power index a (l \ 
lg(M/M ),i.e. 

^m,m q) MM/M & )] oc 



or a Gaussian distribution in 



cxp {-0.5 MM/Mq) - \g(M^/M Q )] 2 /(a^) 2 } (B2) 



r(i) 



and the width of the distribu- 



with the characteristic mass Aft 

( i) 

- . The proportionality factors thereby have to be chosen so 



tion a 

that the mass spectrum $i s (m/m, 



(i+l) 



© 



is continuous at all Af, 



(i) 



and its integral over all \g(M / Mq) becomes unity. The 
choice for the parameters corresponding to different select ed mass 
ranges for disk or bulge lenses following IChabrieJ J2003I) . that is 



adopted for this paper, is shown in Table lBll and the corresponding 
mass function $i g (M/Mg) Q&(M/Mq)) is shown in Fig. lBll 

B2 Mass density 

The view from the observer to the source is reflected by a coor- 
dinate system with a basis (e x , e y ,e z ) where e x points from the 
observer to the source, while e y and e z span a plane perpendic- 
ular to the line-of-sight. For describing properties of the Galaxy, 
however, it is more appropriate to refer to the galactic coordi- 
nates (D^,l,b) which are the spherical coordinates that refer to 
the basis (Ex, Ey,Ez), where Ex points to the Galactic centre 
(l,b) = (0°, 0°), Ey towards the direction of local circular motion 
(l,b) = (90°, 0°), and E z towards Galactic North b = 90°, so 
that X = Dl (cos I cos b, sin I cos b, sin b). The basis (e x ,e y ,e z ) 
arises from (Ex,Ey, Ez) by a rotation TZ(Ez',l) around Ez by 
the angle I and a subsequent rotation 1Z(— lZ(Ez; l)Ey', b) around 
—TZ(Ez; l)Ey by the angle b, i.e. ej = . b)Ej, where 



T(l,b) 




(B3) 



so that vector components transform as x — (T 1 ) (7, b) X = 
T(l,b)X or X = T- 1 {l,b)x = T T (l,b)x. 

The density of matter in the Galaxy is more easily expressed 
in coordinate frames whose origins are at the Galactic centre rather 
than at the position of the Sun. Just by subtracting the correspond- 
ing difference location vector, one finds for such a system with the 
fractional lens distance x = Dl/Ds and the galactic coordinates 
(I, b) of the source star: 

Ro, 



X = X — Ro — x Ds cos I cos b 
Y = Y — x Ds sin I cos b , 
Z = Z = x D s sin b . 

Hence, the cylindrical distance from the Galactic centre is 
R 



= \/x 2 + Y 2 
= Ro \J [cos I - 



x(D s /Ro) cos 6] 2 +sin 2 ^ 
while the spherical distance reads 



(B4) 



(B5) 



v¥ 2 



Y 2 + Z 2 
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Figure B2. Weighted probability density n <fr x (x) of the fractional lens dis- 
tance x = Dl /Ds for two different positions of the source, and lenses 
in the Galactic disk or bulge. The weight factors are given by K dlsk = 
S disk /(S disk + £ bul s c ) and K bul s c = £ bul s c /(£ disk + £ bul s c ), 

so that k& x (x) reflects the mass density, where J [K dlsk< 3> dlsk + 

K bulg C<I) bulgej dx = . e K di S k + K bulgc = ± While K disk ^ 2 x 10 _5 

for the source in the spiral arm, one finds K dlsk = 0.33 and K bulge = 0.67 
for the bulge source. 



= y/Rl + x 2 Dl - 2 xD s R cos I cos b . (B6) 

As chosen by iGrenacher et alj 1 19991) . the mass density of 
the disk is modelled by two double-exponential disks, following 
Bahcall, Soneira & Schmidt 1 1983) and Gilmore, Wvse & Kuiiken 
ll989l) . where~ 

R- iV 



p diak (R,Z) 



0.5 exp 



flthin 



exp 



\Z\ 



Hthin 



. Sthick 

+ — exp 

tl thick 



Ha 



(B7) 



with h = 3.5 kpc being the scale length in the galactic plane, while 
Hthin = 0.3 kpc and //thick = 1.0 kpc are the scale lengths of 
a thin and a thick disk perpendicular to the galactic plane, and the 
corresponding column mass densities are E t hi n = 25 Mq pc -2 
and E thick = 35 Mq pc~ 2 . 

Similar to the di scussion of lHan & Gould ll995blal) and 
IGrenacher et alJ 1 19991) . l et us adopt a mo del of a barred bulge 
based on the COBE data <Dwek et allll995h . which is tilted by an 
angle 9 = 20° with respect to the direction of the Galactic centre, 
so that coordinates along its main axes with origin at the Galactic 
centre are given by 

X' = Xcos6 + Ysm6 y 
Y' = -Xsin6> + Ycos(9, 

Z' = Z. (B8) 

In these coordinates, the mass density of the Galactic bulge can be 
expressed by 



P 



bulge 



(X',Y',Z' 



Po 



bulge 



exp{-s 2 /2} , 



(B9) 



(X>/a) 2 +(Y>/b) 2 +(Z'/c) 



(B10) 



0.62 kpc, and c 

-,10 



0.43 kpc. 

bulge 



with a = 1.58 kpc, b 

A total mass M bu ige = 1-8 X 10 lu Mq implies 
M b uige/(6.577T abc) = 2.1 x 10 9 M kpc~ 3 . 

For a source in the Galactic bulge towards Baade's window 
at (l,b) = (1°,— 3.9°) at a distance Ds = 8.5 kpc as well as 
for a source in the Carina spiral arm at (I, b) = (290.8°, —0.98°) 
and a distance Ds = 6.8 kpc as example for an off-bulge 
target, the weighted probability densities k$ x (x) of the frac- 
tional lens distance x = Dl/Ds are shown in Fig. IB2I The 
weight factors K disk = S disk^ E disk + E bui gc - ) and K bui gc = 

E bulge / ,( E disk + E but ge ) haye been chQsen SQ that Jl j K disk $ disk + 
K bulg C(1) bulgo-| dx = 1 ^disk + K bulgc = L N()t surprisingly! 

the contribution of bulge lenses is negligible for a source in the 
spiral arm, where for the chosen parameters, 

K bui g o = 2 x 10 - 6i 

In contrast, for the bulge source, one finds contributions of com- 
parable order, where K disk = 0.33 and K bulgc = 0.67. While the 
lens mass density for the source in the spiral arm shows a broad 
distribution favouring smaller lens distances, one finds a moderate 
increase with distance for disk lenses and a source in the Galactic 
bulge, while bulge lenses yields significant contributions only for 
x > 0.6. 



B3 Effective transverse velocity 

The effective transverse velocity in a plane at the lens distance 
Dl = x Ds perpendicular to the line-of-sight is given by 



v(x) = i>l — x vs — (1 — x) vq . 



(Bll) 



where «l, i>s, and vo denote the perpendicular velocities of the 
lens, source, or observer, respectively. Let us consider an expec- 
tation value v° = (v), and the source and lens velocities follow 
Gaussian distributions, where isotropic velocity dispersions are as- 
sumed for both the Galactic disk and bulge. While the introduction 
of anisotropies heavily complicates both the discussion and the cal- 
culation, the results are only marginally affected, and the arising 
differences do not exceed those resulting from uncertainties in the 
velocity dispersions themselves. Discussions of anisotropic veloc- 
ity dispersions must not miss the non-diagonal elements of the ve- 
locity dispersion tensor for directions that do not coincide with the 
main axes of the velocity dispersion ellipsoid. 

In this paper, values of cr dlsk = 30 kms -1 for the Galac- 
tic disk and o- bulgc = 100 kms -1 for the Galactic bulge 
have been adopted. Therefore, for bulge sources, the total ve- 
locity dispersion is a(x) = y/l + x 2 a - bulge for bulge lenses or 
u(x) — y/ 1 x 1 (<j bulgc ) 2 + (<r dlak ) 2 for disk lenses, while for disk 
sources, where the lens also resides in the Galactic disk, a(x) = 

VI +X 2 <7 disk 

The probability density (v) of the absolute effective veloc- 
ity therefore takes the form 



$v(v,x) 



2n[a(x)] 2 



271 [a(x)] 2 



exp 



v 2 + [v°(x) 
<X "^ 2[a(x)} 2 ' } ' 



24 M. Dominik 



exp 



v 2 + [v°(x)} 2 \ f vv°(x) \ 



(B12) 



with <(3 being the angle between v and v°(x), v — \v\, v°(x) = 
and Io(r/) denoting the modified Bessel function of the 
first kind to the order zero. 

Hence, with dimensionless £ = v/v c , £ (a;) = v (x)/v c , and 
<r(a;) = a(x)/v c , one finds 

*c(C.*) 

c 



[*(*)]* 



exp 



^ + ic»n l0 (cm\. (B13) 



2[a(x)]2 



[*(*)] 



While the direction of the velocity vector for Bulge objects 
is purely random, disk lenses as well as the Sun perform a sys- 
tematic rotation around the Galactic centre with a velocity v clIC (R) 
depending on the cylindrical distance R. With R given byEq. <B5l . 
the systematic lens motion reads 

, Ro cos I — x Ds cos b 

Z VciTc(R) 



/ 

V-L,y(X 



v °A x ) 



Vclvc(R) 



R 

Ro sin I sin 6 
R 



(B14) 



The rotatio n velocity can be effective l y desc ribed by the model in- 
troduced bv lNavarro. Frenk & White) ll997h . where the mass den- 
sity is given by 



p(x) oc 



1 



r(r + rs) 2 
so that with 



M(r) =4rt / p(r)f '■ df 
o 

andu c irc(r) = [G Af(r-)/r] 1//2 , one finds 

Voire (-R) = Wcirc(-Ro) /nFW X 



(B15) 



(B16) 



(B17) 



where the choices Vd TC (Ro) = 220 km s" 1 for the reference dis- 
tance Ro = 8.5 kpc and rs = 20 kpc yield /Sfw ~ 4.23. 

With respect to the rest frame of the Galaxy, the Sun, located 
at a distance Ro from the Galactic centre, shows a peculiar motion 
V Q = (Vq iX ,Vq jY ,V®,z) = (9,12,7)kms- 1 on top of the 
circular motion of the Galactic disk of Vq circ = (0, v a i IC (Ro), 0) 
with Ucirc(-Ro) = 220 kms -1 . 

One also might consider including the velocity of the Earth 
of v® = 30 kms -1 . While for event time-scales ts 1 yr, this 
velocity is approximately constant (and roughly equivalent to the 
value at the closest angular approach between lens and source), and 
to next order, the acceleration of the Earth's motion can be included 
in the model of the observed light curve by means of additional 
parameters, the full annual modulation affects the light curve for 
longer time-scales and this parallax effect needs to be accounted 
for. In the last case, however, there is no effective Earth's velocity 
that contributes to v. For the calculations in this paper, the velocity 
of the Earth is neglected, so that with T(l, b) from Eq. IB3I . one 
obtains v% =T(l,b)(VQ >cixc + Vq) or 
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Figure B3. Probability density $i g f(lgC) = C m 10 *c(0< wnere 
f = v/vc • The two panels correspond to different positions of the source, 
where the lens is either located in the Galactic disk or bulge at a fractional 
lens distance x = D^/Dg. The bulge population is not considered for the 
source in the spiral arm, because its contribution to the event rate can be 
neglected. 



Vo,!i — — cos / sin b Vq x — 
— sin I sin b [Vq ^ +v c itc(Ra)\ +cos6Vq z . (B18) 

Fig. lB3l shows the distribution of the effective velocity for a 
source in the Carina spiral arm or in the Galactic bulge for either 
bulge or disk lenses, where the same parameters as for the distri- 
bution of the lens distance shown in Fig. IB2I have been adopted. 
In consistence with the latter, 'typical' values of x = 0.85 for the 
source in the Galactic bulge or x = 0.35 for the source in the spiral 
arm have been chosen. The shift towards larger velocities for bulge 
lenses as compared to disk lenses for a bulge source reflects the 
larger velocity dispersion of the bulge, whereas the smaller typical 
velocities for the source in the spiral arm result from the smaller 
velocity dispersion for disk sources and lenses. 



APPENDIX C: STATISTICS OF BINARY ORBITS 

In general, galactic microlensing light curves only depend on the 
components of the orbital separation of a lens binary that are per- 
pendicular to the line-of-sight. Moreover, if the orbital period is 
sufficiently large as compared to the duration of the event, only the 
actual projected orbital separation f rather than the semi-major axis 
a is relevant, where a best-fitting model parameter d — f /rs can 
be determined from the collected data. However, one is interested 
in statistical properties that refer to a rather than to f. For a given 
orbital numerical eccentricity e, the absolute value of the orbital 
separation is given by 



r(t) 



a(l 



1 + e cos tp 



(CI) 
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Figure CI. Probability density 3>x(x) °f the projection factor \ = f /a 
between the semi-major axis a of the orbit and the actual projected sep- 
aration f perpendicular to the line-of-sight. For exoplanets, a distribution 
<J> e = (4/7t) Vl — e 2 has been assumed and the arising results are com- 
pared with the assumption of circular orbits. 



x Q(x-x) 



(C4) 



where the area of the hemisphere (2 n) cancels out against the inte- 
gral over the azimuthal angle (p. 

For a given orbital eccentricity, one obtains the probability 
density of \ as 

oo 



x / , e(x-x)e[x-(i-e)]e[(i + e)- x ] 



71 J VIX 2 -X 2 } [x-(l-e)] [(l + e)-x] 
v 







7t ^/(l + e-x) (1-e + x) 



if 2 



(l + e-x) (1-e-x) J 

for x ^ 1 — e 



= < 



1 /(l + e-x)(l-e + x) 



1 \/- #1 o > 

n y e \ 2 \f ex 



(C5) 







for 1 — e<x<l + e 
for x ^ 1 + e 



so that for the orbital eccentricity e being distributed following the 
probability density <J> e (e), the probability density of \ results as 



$ e (e) $ x (x;e) de, 



(C6) 



where r m i n = a(l — e) and r max = a(l + e) are the mini- 
mal and maximal separations corresponding to the phase angles 
<p = or ip — 7t. With P denoting the orbital period and 
Wcirc = (2n/P)a, one moreover finds the maximal velocity 
«max = fcirc + e )/(l — e) occuring at the minimal sepa- 
ration and the minimal velocity v m - ln = v c i IC — e)/(l + e) 
occuring at the maximal separation. The conservation of angular 
momentum then yields 



, ., 2 dip 2n 2 rr 



(C2) 



Therefore, the probability density of x = r /a, relating the semi- 
major axis a and the actual separation r reads 



P/2 



*x(P;e) = s 



r(t) 



while for circular orbits, one obtains 



*x(x) = 



:G(1- X ). 



(C7) 



Both for circular orbits and elliptical orbits that correspond to 
planetary systems, 3> x (x) is shown in Fig. lCll For the latter case, 
$ e (e) = (4/7t)\/l — e 2 has been chosen in rough agreement with 
radial velocity searches, where p e (e) is approximately constant for 
moderate e, but drops off to zero as e — > 1. 



7ta 2 \/l — e 2 



X 9[x-(l-e)] 9[(l + e)-x] 
71 VIX - (1 - e)] [(1 + e) - x] 



(C3) 



which becomes $ x (x) — 5(x — 1) f° r a circular source, for which 
there is a constant orbital radius r = a. 

An isotropic orientation of the orbit means that the position of 
the companion from the primary at a given phase is uniformly dis- 
tributed on a hemisphere with radius r, so that a probability density 
°f X = r/a for a given \ = r/a reads 



tt/2 



$|(x;x) 



5(p — p cos 0) cos 9 A8 



